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, Abstract 

We consider stochastic dynamics for a spin system with mean field interaction, in which the 
interaction potential is subject to noisy and dissipative stochastic evolution. We show that, in the 
thermodynamic limit and at sufficiently low temperature, the magnetization of the system has a time 
periodic behavior, despite of the fact that no periodic force is applied. 
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1 Introduction 

Dynamics of stochastic systems with many interacting components are often described in terms of the 
interaction energy H(x) — H(xx, X2, ■ ■ ■ , %n)i by stochastic equations of the form 

dX t = -VH(X t )dt + adB t , (1.1) 



where a is the diffusion constant and B is a Brownian Motion. The collective N — > +oo behavior depends 
on the entropy- energy balance: the stochastic term o~dB t in (1.1) is a source of disorder, and it competes 
with the drift —WH(X t ) which tends to "freeze" the system in the states of minimal energy. 

This structure, which admits various versions with discrete time and/or space, has proved to be suc- 
cessful in many applications, starting from physics but then developing in other fields such as Economics, 
Biology and Engineering. 



A substantial modification of the structure above consists in introducing an interaction energy which 
is not a given deterministic function of the state, but it has its own (possibly random) evolution. Models 
of this type have been recently proposed, independently, in various fields of research, but in particular 
£f) • in Biology and Economics. A first example (see e.g. [18]) is provided by the following modification of 

(1.1), that we write componentwise: 



dA7 = d x h(Xl,t)dt + adB l t (1.2) 

N 

d t h(x,t) = -ah(x,t) + Dd 2 x h(x,t) + P^g(X l t ,x). (1.3) 

i=l 

This system has arisen, for instance, in cellular dynamics where cells tend to move toward positions x 
of high concentration h{x, t) of a given chemical (equation (1.2)). This chemical is subject to dissipation 
(a > 0) and diffusion (D > 0). Moreover, while moving, cells release themselves chemicals into the 
medium, thus modifying the concentration h{x,t). It should be noticed that the model is of mean field 
type, in the sense that it is left invariant (in distribution) by permutations of cells. This is one of the 
instances where mean field models, which in physics are often regarded as toy models, play a more 
substantial role. A discrete-time version of this model has been recently studied in [1] . 

Cellular dynamics provides other examples of interacting systems. A very active field of research 
is that of synchronization in synthetic biology (see e.g. [5,12,16]). Concentrations of certain molecules 
within cells vary over time due to chemical reactions, and may lead to periodic behavior. If x denotes a 
vector of concentrations referred to a single cell, its evolution can be described by an ordinary differential 
equation 

x = F{x) (1.4) 
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which exhibits limiting cycles. Given a system of N cells, various interaction mechanisms have been either 
observed or reproduced in experiments. As shown in [5], suitable molecules can be "pumped" into the 
extracellular space, reach through the cell membrane the intracellular space and influence the dynamics 
of x. Denoting by X{ the concentration vector of the i-th cell, by S the extracellular concentration of the 
injected substance and by s, its concentration within cell i, (1.4) is modified as follows 



where s := ^ J^i s i- Random noise can be added to these equations: in all cases the interaction is driven 
by an "external" variable S, which is subject to dissipation (a > 0), and it is of mean-field type. 

A third example we mention has been proposed in [7], it is related to finance and it models defaults 
in a network of firms. For i = 1,2, . . . , N , let <7j(i) € {0, 1} be the indicator of the i-th firm; in other 
words, <7j(t) = means that the i-th firm has defaulted. The dynamics is assigned by giving the 
default rates (probability of default per unit time) A,(t). To model clustering of defaults, one can choose 
Aj = Aj(<ri, . . . , <tjv) to be a deterministic function of <7i, . . . , <jn, decreasing with respect to the natural 
partial order in {0, 1}^ (see e.g. [3]); in this way, the default of one firm permanently increases the rate 
of default of other firms. As observed in [7], the permanence of this effect is unrealistic, and should be 
subject to a sort of dissipation. They propose the following stochastic dynamics for the default rates: 



where \ > are given reference rates, L^{t) := o~i(t), the Bi's are independent Brownian motions, 

and X is some driving exogenous factor, such as a macroeconomic index. If A^ does not depend on i, 
this is again a permutation invariant model, subject to dissipation. 

In terms of rigorous analysis, some results concerning the above models are available. In particular a 
macroscopic equation, which describes the limiting behavior in the limit N — > +oo, has been derived for 
a discrete-time version of (1.2), (1.3), together with sufficient conditions for this macroscopic equation 
to have a globally stable fixed point [1]. For the system in (1.6), a macroscopic equation has also been 
derived, together with a Large Deviation analysis of the N — > +oo limit. 

In comparison with what has been obtained for more traditional models motivated by statistical 
mechanics, one aspect that is missing is the detailed analysis of the phase diagram for some special model. 
This could reveal, for instance, what effects dissipation has in the long-time behavior of the system. The 
aim of this paper is to perform this analysis for a simple model that, however, exhibits some of the main 
features of the models above: the interaction is of mean-field type, and it is subject to dissipation. In 
absence of dissipation, this model reduces to (a slight modification of) the well known Curie- Weiss model, 
whose macroscopic properties are well understood for both the statics and the dynamics. In particular, 
for sufficiently low temperature, the system self-organizes, in the sense that spins tend to align producing 
a macroscopic magnetization that becomes constant in the long-time limit. We will see that in presence 
of dissipation this picture changes: self-organization is still present, but the macroscopic magnetization 
fluctuates periodically driven by a two-dimensional ordinary differential equation similar to that of the 
classical Van der Pol oscillator, in which the large damping regime correspond to low dissipation in our 
model. Despite of the simplicity of the model, some analytic aspects are considerably harder than in the 
standard Curie- Weiss model. 

The emergence of self sustained periodic behavior (i.e. not induced by external periodic forces) has 
been recently studied in several models. In particular, the neural networks in [14, 15] and the active 
rotators in [6] appear to be close in spirit to the simpler model we propose. Precise analogies emerge, 
in particular, with the active rotators model. At macroscopic level, the "inactive" system has a one 
dimensional, compact stable manifold of fixed points. Introducing a small forcing term, a slow, possibly 
periodic motion on the slightly deformed stable manifold is introduced. In our model, in absence of 
dissipation and with a suitable choice of variables, fixed points form a one dimensional, non-compact 
manifold, with an unstable part and two stable branches. A small dissipation introduces a slow motion 
of the stable part of the manifold; as the unstable part is met, the motion is driven to the other stable 
branch. This scheme is then repeated, producing periodic motion. 

We finally mention that in [11], various models for noise- induced behavior are presented; in partucular, 
the model we propose has several analogies with the FitzHugh-Nagumo model [4, 13]. 

In Section 2 we define the model, and derive its macroscopic evolution equation. In Section 3 we 
study the long-time behavior of the macroscopic equation in the case the evolution of the spin-flip rates 



Xi = F(x l )+g(s i ) 

Si = ip(xi) - fisi +7(5 - Si) 

S = -aS + k(s-S), 



(1.5) 



d\i(t) = -a(Xi(t) - X l )dt + (3dL N (t) + a^/X~(i) dB^t) + jX^t) dX(t), 



(1.6) 
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is dissipative but not driven by additional noises. Numerical simulations for a more general case are 
given in Section 4. 



2 Microscopic and macroscopic evolutions 
2.1 The microscopic model 

Let s — (si, S2, ■ ■ ■ , Sn) € {—1,1}^ be a configuration of TV spins. We define, to begin with at an 
informal level, a stochastic process (s(t)) t >o by assigning (besides an initial condition) the spin flip rates. 
For s e { — 1, 1}^ and i = 1,2, ... ,N, define s l to be the configuration obtained from s by flipping Si 
(i.e. s\ = —Si) and with all other spins left unchanged. At a given time t > 0, if s(t) = s, then each 
transition — > — Sj occurs at rate 1 +tanh(sjAj), where Xi(t) is itself a stochastic process, which evolves 
according to the stochastic differential equation 

dAj(t) = -a\i(t)dt + adBi(t)+d(f>(a(t)), (2.1) 

where a, a > 0, 

« w 

^--/Smj^a) (2.2) 

fe=i 

with /3 > 0, and B\,B 2 , ■ . ■ , Bn are independent Brownian motions. A different <fi could be chosen to 
generate other models; we concentrate on this special choice. Note that between two consecutive spin 
flips, the Xi's evolve as independent Ornstein-Uhlenbeck processes; at a spin flip time t, each A^ jumps 
by a quantity of 4>(s(t j) — (f>(s(t—)). 

This construction can be made rigorous. In particular, it defines a Markov process (s(t),X(t)) e 
{—1,1}^ x R N whose infinitesimal generator is given by 

JV 



L N f(s, A) := J2 [(! + tanhM,)) (/(s\ A + lV^s)) - f(s, A)) 



d_ 



-a\i — f{s,X) + —^y2.f(s,X) 



, (2.3) 



where 1 is the TV-dimensional vector with each component equal to 1 and Vi<j)(s) := (j)(s l ) — <f>(s). 

For a clearer picture, consider the case in which a = a = 0. In this case, equation (2.1) is easily 
solved: Xi(t) = c — /3mjv(i), where c = A +/3mAr(0). It follows that the process s(t) is a Markov process 
with infinitesimal generator 

N 

L N f(s) = ]T [1 + tanh( Sl (c - (3m N (s)))} [/(**) - /(*)]. 

i=l 

This is a modification of the standard stochastic dynamics of the Curie- Weiss model, where, in the usual 
version, 1 + tanh(si(c— /3mjv(s))) is replaced by exp[sj(c — /3m;v(s))], where /3 is the inverse temperature 
and c is the magnetic field. As a > 0, in the dynamics of Aj a dissipation, i.e. attraction to zero, is 
added to the jumps; as a > 0, a noise in this dynamics is introduced, in addition to the Poisson-type 
noise driving the spin-flips. 

remark 2.1. For the spin flip rates, that we set equal to 1 + tanh(sjAi), other choices could be made, for 
instance exp[sjAj], that would lead to similar results. The boundedness of our spin flip rates is convenient 
for the proofs, but we believe it is not an essential ingredient. 

2.2 The macroscopic equation 

We now derive the dynamics of the system in the limit as N — >• +oo. In this section we proceed at a 
heuristic level; a formal proof will be sketched in the Appendix. 
We introduce the following empirical measure 



1 N 

^-jvE^W- (2-4) 
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When s = s(t), A = A(t), we write for P N . For an arbitrary function / : {—1, l}xl->l, empirical 
averages will be written in the form 



N 

(p t N j) ■= J fdp? - ^E/^W'^w) • 



The infinitesimal generator (2.3) applied to a function $ of the form 

9(P N ):=4>{{P N ,f)) (2.5) 

reads 

N 



L N c^((P N J)) = ]T [(1 +tanh( Sj A,)) (0«P^",/» - ^((P^,/))) 



9 , 0~ 



(2.6) 



with 
where 



2/3s- 

«i = *fc(i - - *jk*fc . A i = A fe + ^r ( 2 - 7 ) 



TV 

are the spins and the A's after the jump of Sj. 
If we expand the generator (2.6) for large TV we obtain 



C N d>((P N ,f)) « £(1 + tanh( Sj A J ))0'((P JV , /)) 



^(P^, 9a /) + 1(5 ( _ S3 , A3 . ) -^, Aj . ) ,/)) 



+ 0'((P Ar ,/))(-«(P Ar ,Aa A /) + y(p N ,al/)j+o(^-j . (2.8) 

This allows to identify the weak limit P t of the empirical measures , evaluated along the paths of the 
process (2.3), as the (deterministic) solution of the equation 

(P t ,f) - (P , /) = f (P u ,L(P u )f) du , (2.9) 
Jo 

where 

C(P t )f(s,X) = ((l + S tanh(A))(/(- s ,A) - /(a, A)) 

2 \ 

+2f3{P t ,s + tcmh(\))d x f(s,\)-a\d x f(s,\) + ^-d 2 x f(s,\)\ . (2.10) 

We remark that the operator (2.10) can be associated to the nonlinear Markov Process (see [10]) 
{(St, A t )} t>0 , with E t e { — 1, +1} and A t e M, solution of the stochastic equation 

T t — > — T t with intensity 1 + tanh(S t A t ) , 

dA t = (-a A t +2(3(P t ,s + tanh(X)))dt + crdB t , (2.11) 
P t = Iaw(E t ,A t ) . 

Alternatively, writing formally P t = p t (s, A)dA, p t can be identified as the weak solution of the 
nonlinear equation 

d t p t (s,X) = (1 - stanh(A))pt(-s,A) - (1 + stanh(A))p t (s, A) 

2 

+ U — d 2 x p t (s, A) - 2pg(t)d xPt (s, A) + ad x (\ Pt ( S , A)) , (2.12) 
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where 



Pt (s,A)( S + tanh(A))dA 

-co 



The regularizing effect of the second derivative guarantees that p t is indeed smooth in A. 
Notice that defining 



s=±l s=±l 

- and so A) = ^+ u(x \ A) = zH^Hi^ - you get from the Fokker-Planck equation (2.12), 

the equivalent following system of PDEs 

2 

d t v t (\) = ^-d 2 x v t (X) + (aX - 2(ig{t))d x y t {\) + av t (X) , 

d m (\) =^dliit{\) + (a\ - 2/3g(t))d xf H(\) + (a- 2) Mt (A) (2A4) 



- 2tanh(A)f t (A) 



where 



g(t) = (v t , tanh(A)) + (/it, 1) = (i/ t) tanh(A)) + m(i) , (2.15) 

with m(t) the mass of \i at time t or, equivalently, the expected value of S t (the 'average spin' or 
magnetization). Clearly, the knowledge of the pair (^ t ,/i t ) of measure- valued processes, is equivalent to 
that of P t . By definition, v t (X) is the density of the marginal distribution A t , while fi t is the density of 
a signed measure; they must satisfy J R ft (A) dA = 1, J R /it* (A) dA = m(t) G [—1, 1]. 

All the above can be translated into the following rigorous statement, whose proof is given in the 
Appendix. In Theorem 2.1 below, we suppose that the laws of (a N (0) , X N (0)) are P -chaotic. Recall the 
definition of chaoticity. Let 9 be a probability measure on a Polish space X and, for N e N, let Qn be 
a symmetric probability measure on the product space (the law of (a N (0) , X N (0)) is a probability 
measure on ({ — 1, 1} x R) N , assumed to be symmetric). Then (Q^NeN is said to be ^-chaotic if for 
every n € N the joint law of the first n marginals of Oat converges weakly to the product measure ® n 9. 

Theorem 2.1. Assume that the initial conditions (er(0),A(0)) = (a N (0) , X N (0)) for the processes (2.3) 
are Po _cria otic for some probability measure Pq on { — 1, 1} x R. Then the sequence of measure-valued 
random variables (P Ar )iveN converges in distribution as N — > +oo, in the topology of weak convergence 
of probability measures, to the law P on path space of the unique solution of equation (2.11) with initial 
distribution P ; moreover, (Pt)t>o, the measure- valued process of time marginals of P, solves the integral 
equation (2.9). 



3 The case without noise 

In this section we analyze equation (2.9) in the special case of a = 0. We also assume that the initial 
condition is such that all the spins have the same A, i.e. we take Aj(0) = Ao for all j's. These simplifi- 
cations allow a detailed analysis of the long-time behavior of the solution of (2.9). Indeed, writing P in 
the form 

P (s,dA) =5{X-X ) xd Q (s) , (3.1) 
where do is a probability on { — 1, 1}, the corresponding solution of (2.9) maintains the same form: 

P t (*,dA)=*(A-A(t))xd t (s) ) (3.2) 

Plugging (3.2) into (2.9), and setting m(t) := d t (l) — d t (—l) = (d t ,s), we have that the measure (3.2) is 
indeed a solution of the limiting dynamics (2.9), provided that 

A(t) - 20(m(t) + tanh(A(t))) - aX(t) , 
rh(t) = -2(m(t) +tanh(A(t))) ; 

with A(0) = Ao and m(0) = (do, s). 

Notice that, in terms of the pair \xt,vt introduced in (2.13), the solution (3.2) reads 

v t (X) = S(X - A(t)) , jUt(A) - m{t)S{X - A(t)) . (3.4) 



5 



Analysis of the attractors 

First of all, we note that the only fixed point in (3.3) is the origin (0,0) of the phase plane to, A. The 
linearization around this point gives 



m\ _ f-2 -2 \ /to 
Ay I 2/3 2(3 -a) I A 



(3.5) 



and the eigenvalues of the system are 



a 



•'■ = - -± \/| ; - ' - 2j 2 - 2 «- (3-6) 



These eigenvalues have both negative real part for (3 < § + 1, and both positive real part for /3 > ^ + 1. 
Thus, for /3 > § + 1, the local stability of the origin is lost. Much more than local stability can be 
obtained for system (3.3). 

Theorem 3.1. (i) For /3 < f + 1 the origin is a global attractor for (3.3). 

(ii) For (3 > ^ + 1 the system (3.3) has a unique periodic orbit, which attracts all trajectories except 
the fixed point. 

Proof. It is useful to perform a simple change of variable, consisting in replacing to by y := 2(A + (3m). 
In the variables (y, A), system (3.3) becomes 

y = — 2aX, 

(3.7) 

A = y - .9(A) , 

with g(X) = (2 + a)X — 2/3tanh(A). The system (3.7) is of the Lienard type (see, for example, [2, 17]), 
which allows a detailed study of the global stability. 

Case /3 < f + 1. In this case it is easy to show that the function g is strictly increasing, and it is odd. 
Setting 

W(X,y):=^ + ^, (3.8) 

we have 

W = —Xg(X) , (3.9) 

where 



W(X,y) = jW(X(t),y(t)) 



t=o 



and (X(t),y(t)) solves (3.7) with (A(0),j/(0)) = (A, y). Thus W is a global Lyapunov function, which 
implies global stability of the origin. 

Case (3 > § + 1- The odd function g has now two additional, symmetric zeros ±A*, A* > (see figure 1). 
Thus, the system (3.3) satisfies the condition of Theorem 1.1 in [2], which establishes existence and 
uniqueness of a globally stable periodic orbit. 

Here we briefly sketch a proof of the existence, uniqueness and stability of the limit cycle. First of all, 
let us prove that all the trajectories revolve around the origin. Suppose to start in region I (see figure 1): 
y is decreasing, A increasing, so you will eventually hit the nullcline y = g (the origin is the only fixed 
point, and it is repelling). So we are in II: both y and A are decreasing, so either you go in III or you 
'die' at infinity. Let us examine the case y — > — oo; for the slope of the trajectory in this region we have 

dy V —2aX „ 

T7 = - = 7TT ► 3.10 

dA a y-g{X) v^-oo 

which means that we will always end up in III. Ill and IV behave the same, by symmetry. 

Denoting with yo (yi) the intersection of the orbit with the positive (negative) y-axis (which we have 

proved to exist for every trajectory), the zeros of the function 

2 2 

AW(y ) := W(0, Vl ) - W(0, y ) = V -±^ , (3.11) 
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Figure 1: Qualitative behavior of the dynamical system (3.7). The four regions indicate the four different 
directions that the vector field assumes. The blue curve is the graph of g. 




Figure 2: Plot of the phase plane of the system (3.7). One (half) trajectory is sketched (dotted curve), 
starting at (Q,yo) and then back to the y-axis at (0, yi). A (B) denotes the point in which the A- 
component of such trajectory becomes greater (lower) than A*. 
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correspond to periodic orbits (we make reference to figure (2) for notations). 



Name y^ the positive y-axis intersection of the orbit that passes through (A*,0); then, calling t\ the 
time in which (0,yi) is reached, you have 

&W(y* ) = f 1 W(t) dt = - f 1 X(t)g(X(t)) dt > , (3.12) 
Jo Jo 

since </(A) > for A e (0, A*). For y < y$ you get an x-axis intersection smaller than A* -for the 

uniqueness of trajectories. Then you still have AW > 0. 

Now take yo > y^- In this case it is convenient to split AW as follows 

AW(y ) = - f B X(t)g(X(t)) dt - j" X(t)g(X(t)) dt - t X(t)g(X(t)) dt . (3.13) 
Jo Jt A Jt B 

(a) (b) (c) 

The first term is clearly positive, and with a change of variable it can be rewritten as 

(3I4) 

We claim that 

(a) \ monotonically as yo — > +oo . (3.15) 

In order to prove this claim, first of all notice that if one starts at (0, 2/0 = 2/0 + e ) f° r e positive 
and arbitrarily small, then y(A) > y(A) for every A between zero and the A-axis intersection with the 
trajectory. This implies that (a) decreases monotonically with yo. Moreover, the slope of the trajectories 

dy _ 2aX 
dA _ ,g(A) - y 

is bounded as long as we are away from the nullcline y = g(X) and on a compact interval in A. So, given 
an arbitrary positive number M you can always find an initial condition (0,yo) an d a A > A* such that 
y(X) > M as long as A e (0, A). This proves the claim. 

Analogously, one can prove that (c) is positive and monotonically decreasing to zero as yo — > +oo. 
Now let us deal with the second term in (3.13). It can be rewritten as 

(&) = - g(X(y))dy , (3.16) 

which is negative. We claim that (b) \ — oo monotonically as yo — > +oo. In order to show this, it is 
convenient to split (b) as follows 

rtA+S rt B —S rt B 

(b)= / Wdt+ / Wdt+ / Wdt, (3.17) 

Jt A Jt A +S JtB— 6 

with (5 arbitrarily small (and positive). The first and third terms in the right hand side of (3.17) remain 
negative and finite in the limit yo — > +oo. For the second one, namely 

rv(tA+S) 

- g(X(y))dy, (3.18) 

Jy(t B -S) 

the integrand is bounded away from zero (which is not true for (3.16)) and y(tA+S) /* +oo monotonically 
when yo — > +oo, which comes as a consequence of the proof of (3.15). This is enough to prove the claim. 
In the end, we have that 

• AW(yo) is positive when y < y^ 

• AW(yo) \ — oo monotonically as y — > +oo 

which prove the existence and uniqueness of the periodic orbit. In order to prove stability, it is enough 
to say that AW > when y < y^ and AW < when y > Vq, where y^ denotes the positive y-axis 
intersection of the periodic orbit. □ 

Summarizing, the simple model presented in this section is obtained from a standard Curie- Weiss- 
type model by introducing a dissipation on the spin-flip intensity. This dissipation does not destroy 
self-organization of the spins. However, the nonzero magnetization produced by the self-organization 
does not converge to a constant value, as in Glauber dynamics for ferromagnets, but rather oscillates 
periodically. 
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4 Numerical results and conclusions 



Here we would like to present and briefly discuss some numerical results regarding the system of partial 
differential equations (2.14), that is the general case with the extra noise in the dynamics of the intensities 
Aj's (2.1), whose qualitative analysis is intended to be subject of future work. The qualitative idea that 
we would like to stress (supported by the numerical tests) is the following: for small enough values of the 
diffusion parameter a the system qualitatively behaves just as in the no-diffusion case we have analyzed 
in the previous section, i.e. the extra noise gives just small perturbations around the periodic orbit (for 
the super-critical case /3 > a/2 + 1) or around the totally disordered configuration (when (3 < a/2 + 1); 
for cr large enough, instead, the diffusion term dominates and both the super-critical and sub-critical 
cases evolve to gaussian behaviors around the totally disordered configuration (m, A) = 0. 

Figure 3 refers to the sub-critical case, with parameters a = 3, (3 = 1. There it is shown a comparison 
between two phase plots: one is the phase plot of to — (fj,, 1) versus the expected value of A 

(A) t := / xvt(dx) , 

and with a = 0.1, the other one is the phase plot of the corresponding a = case, i.e. it is calculated by 
solving system (3.3). Here it is clear that the trigger of a small (but non-zero) value of extra noise does 
not spoil the qualitative behavior of the system. Moreover, from equation (2.11), it is easy to see that 
the variance of ft (A) 

<(A-(A) t ) 2 )*=Var(A t ) 

satisfies the equation 

-|var(A t ) - -2aVar(A t ) + a 2 , 

so that 

Var(A t ) - e- 2at Var(A ) + ^ (l - e" 2 " 4 ) . (4.1) 
In all simulations the initial condition for v is (approximately) a delta function at Ao = 3, so Var(A ) ~ 0. 

2 

For large t, the variance approaches the value ^j. 

An analogous situation is found in the super-critical case, as can be seen in figure 4 One can see that 
the mean of the densities keeps oscillating in time and that the periodic behavior is indeed preserved, 
even if slightly modified, in the small noise case. By (4.1), the variance of v remains bounded and 
actually quite small during the evolution, thus showing that turning on a small noise does not spoil the 
qualitative behavior of the noiseless system. 

Figures 5 and 6, are the analogues of the preceding ones, but when a = 10, i.e. with a large 
contribution of the extra-noise in the dynamics of the intensities (2.1): it is clear that the diffusion 
is here predominant, pushing the variance up to a gaussian behavior around the totally disordered 
configuration (to, A) = (0,0). The more interesting result is perhaps the one shown in figure 6: the 
periodic behavior of the no-diffusion case is completely lost and the system rapidly evolves to the origin 
of the phase plane. 
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phase plot for m and the average value of X phase plot for the case a=0 




m m 

Figure 3: Phase plot of the expected value of A versus m (left) and phase plot of the no-diffusion case 
(namely of the dynamical system (3.3)) (right). The parameters are a = 3, /3 = 1, a = 0.1, An = 3, 
m(0) = 0. 



phase plot for m and the average value of X phase plot for the case a=0 




m m 

Figure 4: Phase plot of the expected value of A versus m (left) and phase plot of the no-diffusion case 
(right). The parameters are a = 3, (3 — 3, a = 0.1, An = 3 and m(0) = 0. 
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phase plot for m and the average value of X phase plot for the case a=0 




Figure 5: Phase plot of the expected value of A versus m (left) and phase plot of the no-diffusion case 
(right). The parameters are a — 3, (3 — 1, a = 10, An = 3 and m(0) = 0. 



phase plot for m and the average value ofl phase plot for the case a=0 




m m 

Figure 6: Phase plot of the expected value of A versus m (left) and phase plot of the no-diffusion case 
(right). The parameters are a — 3, /3 — 3, a = 10 and An = 3, m(0) = 0. 
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A Proof of Theorem 2.1 



For X a Polish space, denote by T(X) the space of probability measures on the Borel sets of X; equip 
5 > (X) with the topology of weak convergence, which makes it a Polish space, too. As in Section 2, let 
P N be the empirical measure of the iV-particle system. We may assume that the { — 1,1} x R- valued 
processes (sj,Xj) have cadlag trajectories (i.e., trajectories that are right-continuous with limits from 
the left). Consequently, P N is a probability measure on the Borel sets of D := Z?([0,oo),{— 1,1} x R), 
the space of { — 1, 1} x R- valued cadlag functions equipped with the Skorohod topology. 

The strategy of proof, here, is to represent both the microscopic and the macroscopic model as 
solutions of certain stochastic differential equations in order to apply results by [9] on propagation of 
chaos, which implies convergence of empirical measures. 

Let r\ be Lebesgue measure restricted to the Borel sets on the interval (0, 2). Let ((0, 5F, P), (5Ft)t>o) 
be a filtered probability space satisfying the usual hypotheses rich enough to carry an independent family 
(Bi,J^i) ie ^ of one-dimensional (5F t )-Brownian motions Bi and stationary (9 r t )-Poisson random measures 
Kj with characteristic measure r\. For N <G N, consider the system of Ito-Skorohod equations 

fl N f 

dAf (t) = -aAf (t-) dt + adB.it) £ ]T / q ((*?(*-), Af (t-)), u) X k (du, dt), 

N k =l J (0.2) (A1) 

ds?(t)= [ g((*f r (t-),Af r (*-)).«) ^i(d«,dt), ie{l,...,N}, 

J (0,2) 

where 

q((s, A),m) := -2/i(s) • l(o,i+h( s )tanh(A))H, (s,A)eR 2 , ue(0,2), 

and /i(s) := (—1) V (s A 1), s € R. By Theorem 1.2 in [9], existence and uniqueness of solutions hold in 
the strong sense for the system of equations (A.l) since its coefficients are globally Lipschitz continuous; 
the jump coefficient, in particular, satisfies the L 1 Lipschitz assumption of the theorem. Clearly, if 
s G { — 1, 1}, then h(s) = s, h(s) tanh(A) = tanh(s • A), and 

/ q((s, A), u)»7(du) = -2 (s + tanh(A)) . 

J (0,2) 

Thanks to the choice of the jump heights, if (3^(0)^(0)) is such that sf (0) e {-1, 1}, then sf (t) G 
{ — 1,1} for all t > 0. Let us fix a sequence of initial conditions (s N (Q), \ n (0))n£N such that s N (0) £ 
{-1,1} W for all N e N and (P o (s N (0), X N {0))- 1 ) NeN is ^-chaotic for some fj, e ?({-l, 1} x R). The 
solution process (s N , X N ) is then a {—1, 1}^ x R^-valued Markov process. Comparing its infinitesimal 
generator (cf. equation (1.2) in [9]) with equation (2.3) above shows that (s N , X N ) is a realization of the 
A^-particle microscopic model. To prove Theorem 2.1 we thus have to prove convergence of the empirical 
measures P N :— J2iLi $(s N .\ N ) associated with the solutions of (A.l). 
Define a function b: T(R 2 ) -> R by 

6(i/):=2/3- / (ft(s)+tanh(A))i/(ds,dA). 
Jr 2 

Notice that b is Lipschitz continuous with respect to the bounded Lipschitz (or Dudley) metric on J'(R 2 ) 
as well as with respect to the Wasserstcin-1 (or Lipschitz) metric on Ti(R 2 ), the space of probability 
measures with finite first moments. By Theorem 2.1 in [9], existence and uniqueness of solutions hold in 
the strong sense for the McKean-Vlasov Ito-Skorohod equation 

dA(i) = -aA(t)dt + <rdB 1 (t) + b(P t )dt, 

dS(t)=/ g((S(t-),A(t)),«)Ni(d«,dt), (A.2) 

J (0,2) 

P t =Iaw(E(t),A(t)). 

Assume that P = A* = law(S(0), A(0)). Set P := law(S, A) and observe that P E V{D). Comparison of 
infinitesimal generators yields that the solution (E, A) of (A.2) is a realization of the nonlinear Markov 
process given by equation (2.11). Moreover, the ?({— 1, 1} x R)-valued process (Pt)t>o coincides with 
the solution of equation (2.9) with initial condition P - 
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For N e N, consider the system of Ito-Skorohod equations 



dAf (t) = -aAf (t) dt + adBi{t) + b(P t N ) dt, 

(A.3) 

'(0,2) 



dsf(t)=[ g ((*f r (t-),Af(t)),«)N i (d«,d*), ie{i,...,iV}, 

J(0,2) 



where P/^ := ^(A N (t) s N (t)) ^ s em Pi r i ca l measure of the solution at time t. Notice that Xf (t) = 

Af (t— ) by continuity of trajectories and that all processes are stochastically continuous. Again by 
Theorem 1.2 in [9], existence and uniqueness of solutions hold in the strong sense for the system of 
equations (A.3). If the initial condition (s N (0) , X N (0)) for (A.3) is such that sf (0) £ {-1,1}, then 
sf (t) e {-1,1} for all t > 0. Fix the initial condition at (s N (0), X N (0)) := (s N (0) , X N (0)) . Since s N (0) 
takes values in { — 1, 1} N and by the continuity of Lebesgue integrals, the system of equations (A.3) can 
be rewritten as 



8 N r 

dAf (t) = -aAf (i-) dt + a dB t (t) £ £ / g ((*?(*-). A? (t-)), «) tfdu) dt, 

(A4) 

dsf(t)=/ 9 ((sf(t-),Af(t-)),u)N;(d«,dt), ie{l,...,JV}. 

i(0,2) 



Set P w = J^ili ^(s N ,A N )- By Theorem 4.1 in [9], the sequence (law(s Ar , A Ar ))Ar g N is P-chaotic. This 
implies, by the Tanaka-Sznitman theorem (for instance, Theorem 3.2 in [8]), that the sequence (P N )N&i 
of y(P)-valued random variables converges in distribution to the probability measure P. In order to 
establish convergence of (P n )n&n to P, it is therefore enough to show that 



d bL (law(P Ar ),law(P w )) N -±?0 



where dbL is the bounded Lipschitz metric on 7 I (7(D)). By definition of dbL and since both P N and P N 
are empirical measures for processes defined on the same stochastic basis, we have 



1 N 

d bL (law(P^), law(P^)) < E [d bL (P N , P N )] < jj £ E ' X ^ *f ))] 



where df,z, is the bounded Lipschitz metric on 7(D) and dg/co the Skorohod metric on D. For i £ N, let 
Ni be the compensated Poisson random measure associated with Nj, that is, N,(du, dt) = Nj(du, dt) — 
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Tj(du)dt. Then for T > 0, % e {1, . . . , N}, N e N, 



E 



sup |Af(t)-Af(i)| 

t€[0,T] 



< qE 



/■ T |Af(t-)-Af(i-)| 
Jo 



dt 



E 



sup 

te[o,T] 



JV 

E. V 
/v 2^ 



g ((*jS r (r-),A£ r (r-)) ) «)3Sf fc (d«,dr) 



fe=l 



sup 

te[o,T] 



J(0,2) 



k=l J0 J (0,2) 

« ((^(r-), Af (r-)), «) - g ((sf(r-), Af (r-)), «)) „(d«) dr 



<a / T E[|Af(t-)-Af(i-)l]df + 4E A £ f [ q((a?(t-),\Z(t-)),u)X k (du,dt) 

JO JV . , JO J (0,2) 



k=1 J0 J (0,2) 



1/2 



iV 



fe=i 



/ T / |? ((*"(*"), (*-))> u) q ((«?(*-). Af (t-)), u) I d U dt 

JO J(0,2) 



< 



a / T E[| 

Jo 



Af(t-)-Af(t-)|]dt+^E 



( g (( S f(t-),Af(t-)), U )) 2 d« dt 



n 1/2 



+ 



< a 



^ f] / T E [|( S f (t-) - sf (HI + |A^(t-) - Af (t-)|] dt 
fc=i Jo 

/ T E[|Af(H-Af(t-)|]dt + 
Jo 



8/3 V2T 
~7W~ 



R q N l-T 

+ ^Ey o E[|( s f(H-^(HI + |Af(t-)-Af(HI]di. 



Since 
E 



sup | S f (t) - S "f (t)| 

t€[0,T] 



< E 



/O J(0,2) 



< 



9 ((*?(*-), Af (t-)), u) - q ((sf (t-), Af (t-)), «) I d U dt 

6 / T E [|( S f (t-) - sf (t-)| + |Af (t-) - Af (t-)|] dt, 
Jo 



it follows that 
1 N 



i=l 



< 



sup (| s f(t)- s -f(t)| + |Af(t)-Af(t)|) 
te[o,T] 

JV t 



8/3V2T ^ a + 6 + 6/3 



/N 

V U 2 — 1 



iV J 1 

We 0(^(*-) - *T (HI + l A ^(H - *T(HI] d * 
<=i - 70 



T i JV 



sup (|( s f( r )- s -f( r )| + |Af(r)-Af(r)|) 

r£[0,t] 



dt. 



An application of Gronwall's lemma yields, for every T > 0, 



Ive 



sup (| S f(t)- S -f(t)| + |Af(t)-Af(t)|) 
te[o,T] 



which implies the desired convergence. 
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